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X> ' ABSTRACT 

CD ' 

Context. To investigate the properties of the 21-cm radio-lines of galactic neutral hydrogen, the profiles of "The 
' Leiden/ Argentine/Bonn (LAB) Survey of Galactic Hi" are decomposed into Gaussian components. 

Aims. The width distribution of the obtained components is analysed and compared with similar studies by other authors. 
, Methods. The study is based on an automatic profile decomposition algorithm. As the Gaussians obtained for the complex H i profiles 

■ near the galactic plane cannot be directly interpreted in terms of the properties of gas clouds, we mainly study the selected simplest 
' profiles in a limited velocity range. The selection criteria are described and their influence on the results is discussed. 

Results. Considering only the simplest Hi profiles, we demonstrate that for Gaussians with relatively small LSR velocities 
y—{ , (-9 < Vc < 4 km s~') it is possible to distinguish three or four groups of preferred line-widths. The mean widths of these groups are 

. FWHM = 3.9 ± 0.6, 11.8 ± 0.5, 24.1 ± 0.6, and 42 ± 5 kms~'. Verschuur previously proposed similar line-width regimes, but with 

somewhat larger widths. He used a human-assisted decomposition for a nearly 50 times smaller database and we discuss systematic 
| differences in analysis and results. 

Conclusions. The line- widths of about 3.9 and 24.1 kms~' are well understood in the framework of traditional models of the two- 
phase interstellar medium. The components with the widths around 1 1.8 kms~' indicate that a considerable fraction (up to about 40%) 
i-Q ' of the H i gas is thermally unstable. The reality and the origin of the broad lines with the widths of about 42 km s _I is more obscure. 

These, however, contain only about 4% of the total observed column densities. 

O ' Key words. ISM: atoms - solar neighbourhood - Radio lines: ISM 
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1. Introduction be found at all. However, when we decompose a large number 

of profiles by exactly the same algorithm, we may hope that 
The tradition to analyse the 21-cm H i line profiles in terms of for similar proffles the algorithm will prefer similar decompo- 
Gaussian components has its roots in the finding that the H i ab- sitions This means that usually some specific f ea t ur es in the ob- 
£ . sorption features are often well represented by Gaussian func- seryed profiles are repreS ented by some specific set of Gaussians, 
tions in optical depth r(v). This is understandable when the line which can be found from the overall data . set more easily than 
profile is dominated by the thermal motions or nonthermal mo- un . parame trized spectral features. The latter is one of our main 
tions consisting of a large number of turbulent elements. The interests in Gaussian decomposition and from this point of view 
decomposition of the emission profiles is physically meaning- all obtaine d Gaussians are real, as they describe some features in 
ful only when r(v) « 1. Nevertheless, both the opacity and obserV ed profiles. If and how these Gaussians are related to the 
brightness profiles of most sources are easily decomposed into physical prope rties of the galactic ISM is a question, which must 
Gaussains, so whether or not this model is physically correct, it be answere d separately for every selected subset of the compo- 
works empirically and is convenient. Moreover, even if such a nen t s 
decomposition is not a physically valid description of the galac- 
tic ISM, it is a valid mathematical description of the correspond- 
ing 21-cm Hi profile. Proceeding from these considerations, in the first paper of 
Mathematically the most questionable point of the Gaussian this series (Haud 2000 hereafter Paper I) we described a new 
decomposition is that often it is not unique. This nonuniqueness fully automated computer programme for the decomposition of 
means that several quite different solutions may approximate the large 21-cm H i line surveys into Gaussian components. We used 
observed profile almost equally well, and the decomposition pro- the Leiden/Dwingeloo Survey (LDS) of galactic neutral hydro- 
vides no satisfactory means for choosing between these solu- gen (Hartmann |1994| l as test data for decomposition. After our 
tions, while others, equally good or even better ones, may not initial paper a new revisio n of th e LDS data (LDS2) has been 

published (Kalberla et al. [2005 ). Recently a similar Southern 

Send offprint requests to: U. Haud sky high sensitivity Hi survey at 6 < -25° was completed in 

* Founded by merging of the Sternwarte, Radioastronomisches the Instituto Argentin o de R adioastronomia (IARS) and pub- 

Institut and Institut fiir Astrophysik und Extraterrestrische Forschung lished by Bajaja et al. ( 120051) . Both surveys have been combined 

der Universitat Bonn to form the Leiden/ Argentine/Bonn (LAB) full sky database 
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(Kalberla et al 2005 ). Using our Gaussian decomposition pro- 
gramme, we have decomposed all profiles in these surveys. 

The specifications for the LDS2 and the IARS closely match 
each other, but all the data reduction and calibration procedures 
for combining them into the LAB were carried out entirely in- 
dependently (Kalberla et al. 2005). Proceeding from this, also 
the Gaussian decomposition was carried out separately for the 
LDS2 and the IARS. For the published version of the LDS2, in 
the cases of repeated observations at the same sky positions, the 
final profiles were selected on the basis of the best agreement of 
their preliminary Gaussian decompositions with the decomposi- 
tions of the neighbouring profiles (Kalberla et al. 120051 ) and for 
our final decomposition we used these preselected profiles. In 
the case of the IARS, we used for the decomposition the original 
1008-channel data of all the observed profiles and for repeated 
observations we applied before the final decomposition the same 
selection criteria as used for the LDS2. 

In both surveys the observations were made on a galac- 
tic coordinate grid with observed points spaced by (A/, Ab) = 
(0°5/ cos/?, 0°.5). The velocity resolution was 1.03 kms~' in the 
case of the LDS2 and 1.27 kms~' for the IARS. For both surveys 
the authors declare the final rms noise of the data to be about 
0.07 K, but our estimates have given somewhat larger values. 
For the LDS we have obtained an estimate of 0.09 K (Paper I) 
and the same value is given also by Westphalen ( 1 19971 ). With 
the IARS the situation seems to be more complicated (Haud & 
Kalberla 2006, hereafter Paper II). When for the signal-free re- 
gions of the profiles the estimate given by the authors of the sur- 
vey seems to be correct, in the regions containing the emission 
signal, we estimated higher noise levels than it may be expected 
on the basis of the signal-free regions and the radiometer equa- 
tion. From the regions containing the line emission we estimated 
the noise level for zero intensity line emission to be about 0.08 K 
(Fig. 12. of Paper II). 

When performing the Gaussian decomposition, the main at- 
tention was turned on keeping the number of resulting Gaussians 
as small as possible (Paper I). For every profile in the first de- 
composition stage new Gaussians were added until the rms of 
the residuals became less or equal to the noise level of the signal- 
free regions of the profile. The dependence of the noise level 
on the signal strength was taken into account according to the 
radiometer equation. After preliminary decomposition of each 
profile special analysis was applied to find additional possibil- 
ities for removing some less important components from the 
final decomposition without increasing the residuals above the 
noise level. To reduce the ambiguities in choosing between dif- 
ferent, but nearly equally acceptable solutions, we have used 
the assumption that in survey observations the profiles from 
neighbouring sky positions must share some common properties 
and therefore their decompositions must be as similar as possi- 
ble. The described procedure gave us 1 064 808 Gaussians per 
138 830 profiles for the LDS2 and 444573 Gaussians per 50980 
profiles in the case of the IARS. 

In Paper II we analysed the distributions of the parameters 
of the obtained components. We discussed both the original LDS 
data and the newer LAB. The main attention was paid to the sep- 
aration of the Gaussians describing different artefacts of the ob- 
servations (interferences), reduction (baseline problems) and the 
decomposition (separation of the signal from the noise) process. 
In the present paper, we continue the analysis of the distribution 
of the Gaussian parameters, but in this case the main focus is on 
the components, most likely corresponding to the galactic Hi. 
We start with the general distribution of all Gaussians and con- 



centrate on the distribution of the widths of weak components 
with near zero velocities. 

Earlier, a similar distribution has extensively been studied 
by Verschuur and co-authors (fl989l [19941 [19991 [20021 120051 . 
Their results are qualitatively very similar to ours, but we must 
note some differences in the quantitative values. By trying to 
follow the approach used by Verschuur, we point out the main 
differences found and their possible influence on the final results. 
In the concluding sections of the paper, we discuss the relation 
of our results to the properties of the galactic H i. 



2. The general distribution of Gaussian widths 

In Paper II we discussed mainly the frequency distribution of the 
Gaussians in their height vs. width plane (Fig. [TJ. We defined 
the height of a Gaussian by the value of the central brightness 
temperature 7bo > from the standard Gaussian formula 

_(V_v c ) 2 

T b = T w e^V, (1) 

where 7\, is the brightness temperature and V is the velocity of 
the gas relative to the Local Standard of Rest. Vc is the velocity 
corresponding to the centre of the Gaussian. We characterise the 
widths of the components by their full width at the level of half 
maximum (FWHM), which is related to the velocity dispersion 
cr v obtained from our decomposition programme by a simple 
scaling relation FWHM = V8 In 2cr v . 

In Paper II we have demonstrated that not all Gaussians ob- 
tained in our decomposition and represented in Fig. Q] could be 
considered as corresponding to the real H i emission. A consid- 
erable number of the obtained components are due to different 
observational, reductional and decompositional problems. Very 
narrow Gaussians actually represent the strongest random noise 
peaks misinterpreted by the decomposition programme as pos- 
sible signal peaks and the radio-interferences not found or prop- 
erly removed during the reduction of the observed profiles. Many 
of somewhat wider weak Gaussians seem to be caused by the in- 
creased uncertainties in bandpass removal near the profile edges. 
In the direction of even wider Gaussians, at least some weak 
components may be due to the problems in the determination 
of the profile baselines, but some of these Gaussians may also 
arise from the emission of the high velocity dispersion halo gas 
(HVDHG) reported by Kalberla et al. ([19981 ). 

In Fig. [T] such spurious components (in the lower left-hand 
part of the figure) are separated by two thick solid lines from 
the others (in the upper right-hand part of the figure) more likely 
describing the real galactic Hi emission. However, from Fig. Q] 
we can also see that the parameters of the components, most 
likely corresponding to the galactic H i, are not distributed ran- 
domly, but exhibit concentration around two more or less dis- 
tinct values. We can clearly see such a concentration around the 
values of Tbo ~ 8 K, FWHM « 4 kms -1 and a weaker one 
at r b() * 2.5 K, FWHM « 22 kms" 1 . This findin g is in gen- 
eral agreement with the earlier results by Mebold (1972), who 
decomposed nearly 1700 21 -cm Hi emission line profiles into 
about 2400 Gaussian components and found that their width dis- 
tribution has two distinct maxima at FWHM » 7 km s -1 , mostly 
populated by components with 10 ^ Tbo 55 50 K (narrow com- 
ponents), and 24 kms -1 for Gaussians with 3 ^ 7bo 5 10 K 
(shallow components). 

Mebold (1972) used the profiles with lower spectral reso- 
lution than that of the LAB. Moreover, his profiles were from 
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Fig. 1. Frequency distribution of the parameter values in 
(lg(7bo), lg(FWHM)) plane for all Gaussians obtained in the de- 
composition of the LAB. Isodensity lines are drawn in the scale 
of \g(N +1) with the interval of 0.2. The two main maxima are 
labelled with their height values. The green thick solid lines 1 
and 2 represent the selection criteria discussed in Paper II. 

relatively low galactic latitudes and the stronger emission com- 
ponents become wider and more non-Gaussian due to the satu- 
ration and self-absorption. Considering this, we have a surpris- 
ingly good match of the maxima of at first glance physically 
meaningless (lumping into one high-latitude Hi and heavily- 
saturated low-latitude galactic disk profiles, HVC's and IVC's 
and even the observational, reductional and decompositional 
artefacts) distribution in Fig. Q] with the ones described by 
Mebold. In his paper Mebold concluded: "... it is most likely 
that the narrow components are emitted by cold (T ^ 70 K) Hi 
clouds and the shallow components by a hot (750 S T S 9200 K) 
and optically thin Hi intercloud medium". However, in this 
way the general coincidence of the distribution properties of the 
Gaussians described above, may indicate that even in our mathe- 
matical representation of the profiles with the Gaussians at least 
some components may preserve physical information about the 
properties of the underlying galactic H i. In the following we will 
try to separate from the overall set of Gaussians the ones most 
directly related to the structure of the local galactic H i. 

3. The central velocity vs. width distribution 

It is well known that due to the concentration of the gas in 
the Galaxy into a thin disk and the differential rotation of this 
disk, the parameters of the observed H i profiles depend rather 
strongly on the direction of the observations. The area under the 
profile shows the 1 / sin b variation with the galactic latitude. The 
centre of gravity of the line has the sin 21 variation with longi- 
tude and cos 2 bj sin b variation with latitude. The line-width has 
a component which changes like sin 2 21 with longitude and like 
cos 4 b/ sin 2 b with latitude. A considerable part of these changes 
in line-shapes are represented in the Gaussian decomposition 
by decomposing the more complicated profiles near the galactic 
plane with the larger numbers of Gaussians, but most likely to 
a certain extent these dependences will also affect the properties 
of the Gaussians themselves. 

From the relations discussed above, it follows that with the 
increasing mean velocity of the line also the line-width must in- 
crease. For our Gaussians this is illustrated in Fig. [2] separately 



for profiles at latitudes \b\ > 40° (upper panel) and \b\ < 40° 
(lower panel). To construct both panels of this figure, we sorted 
in corresponding latitude intervals all Gaussians with central ve- 
locities in the decomposition range (-460 < Vc < 396 kms~' 
in the LDS2 and -437 < V c < 451 kms~' in the IARS; for 
details see Papers I and II) in the ascending order of their cen- 
tral velocities and then grouped the sequence into 129 bins of 
an equal number of Gaussians, rejecting some Gaussians with 
the most extreme velocities. We binned the line-widths in equal 
steps of lg(o>) of 0.025. The isolines in Fig. [2] give the num- 
bers of Gaussians in each of such two-dimensional parameter 
interval. The latitude limit of \b\ = 40° was chosen rather ar- 
bitrarily. For testing purposes, we constructed similar plots also 
by taking the limit at \b\ - 30°. In this case, the total number 
of Gaussians was increased in the upper panel and decreased in 
the lower panel, but the general structure of the distributions re- 
mained the same. Therefore, we may expect that the exact value 
of the limiting latitude is not very critical until it remains in some 
reasonable range. 

The first thing we can see from Fig. [2] is that the typical 
line-widths of the density enhancements in the distribution of 
Gaussian parameters are not independent of the velocities of the 
Gaussians, but are larger for larger absolute values of the central 
velocities of the Gaussians. Moreover, this dependence seems to 
be stronger for lower latitudes than for higher ones, as expected 
from the relations discussed above. We can also see that at high 
velocities Figs.|2ji (\b\ > 40°) andEfc (\b\ < 40°) are rather similar 
(when we take into account the differences in the velocity scales 
of these plots). This is most likely caused by the fact that even 
near the galactic plane the amounts of the very high velocity gas 
remain modest, the corresponding Gaussians are relatively well 
separated from the bulk of the disk gas and their parameters are 
rather well defined. 

At lower velocities, the differences are more remarkable. As 
discussed in detail in Paper II, there is much more gas near the 
galactic plane than at greater heights and therefore at low lat- 
itudes the effects of velocity crowding, blending and saturation 
become more severe and the profiles become more complex than 
at higher latitudes. As a result, the lower panel of Fig.|2]is dom- 
inated by the narrowest Gaussians and at larger widths we can- 
not distinguish other clear frequency enhancements. Most likely 
such a picture is caused by the fact that in this region the ob- 
tained Gaussians are only mathematically representing the pro- 
file and do not have any direct relation to the properties of the 
galactic H i. However, in the distribution of Gaussian widths at 
high galactic latitudes we can easily distinguish three frequency 
enhancements. Those with the smallest and largest widths corre- 
spond more or less to the enhancements visible also in the gen- 
eral distribution given in Fig.Q] but the concentration at interme- 
diate widths is new, not unambiguously visible in the previous 
figure. 

When the line-widths of the density enhancements in Figs. [2] 
are in general dependent of the central velocity of the corre- 
sponding Gaussians, they are nearly velocity independent in 
the region of -9 ^ Vc ~ 4 kms -1 (marked in Fig. [2^ by 
two solid vertical lines). Therefore, in this velocity interval it 
would be also meaningful to estimate the mean line-widths of 
the three visible line-width groups of Gaussians. To obtain these 
estimates, we first found the distribution of the widths of all 
Gaussians with \b\ > 40°, -9 ^ V c 5 4 km s _1 and modelled this 
distribution as a sum of three lognormal distributions defined by 

C (x < 0) 
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Fig. 2. The distribution of Gaussian parame- 
ters in the central velocity - Gaussian width 
plane for two different latitude ranges: upper 
(a) panel for \b\ > 40° and the lower (b) 
one for \b\ < 40°. The contour lines on both 
panels are drawn at the levels of 1, 2, 5, 10, 
15, 30, 45, 60, 75, 85, 100, 115, 125, 140, 
160, 180, 200, 225, 250 and 300 Gaussians 
per counting bin. The ticks on the x-axis are 
drawn so that the number of Gaussians between 
every pair of neighbouring ticks is approxi- 
mately the same. The precise values of the ve- 
locity corresponding to ticks of the upper panel 
are: -450, -215, -120, -80, -60, -52, -47, 
-41, -36, -30.8, -26.7, -22.5, -19.1, -16.5, 
-14.4, -12.7, -11.3, -10.2, -9.4, -8.7, -8.1, 
-7.6, -7.1, -6.6, -6.2, -5.8, -5.4, -5.0, -4.6, 
-4.2, -3.9, -3.5, -3.16, -2.82, -2.47, -2.13, 
-1.77, -1.39, -0.95, -0.50, 0.00, 0.59, 1.33, 
2.21, 3.2, 4.6, 6.3, 8.6, 12.4, 21.6, 37, 70, 118, 
207 and 440. Those for the lower panel are: 
-450, -138, -107, -89, -74, -63, -54, -48, 
-41, -35, -30.3, -25.7, -22.0, -18.7, -15.8, 
-13.3, -11.3, -9.5, -8.0, -6.7, -5.5, -4.5, 
-3.6, -2.75, -1.97, -1.26, -0.59, 0.00, 0.54, 
1.08, 1.60, 2.10, 2.63, 3.2, 3.8, 4.4, 5.1, 5.8, 6.6, 
7.5, 8.5, 9.7, 11.1, 12.6, 14.5, 16.9, 19.9, 23.6, 
28.5, 35, 43, 52, 62, 73, 88, 107, 156 and 440. 
In panel a, the vertical solid lines enclose the 
velocity region -9 < Vc < 4 kms~' of weak 
dependence of the line-widths on the central 
velocity of the components and the dashed lines 
correspond to the region -31 < Vc ^ 18kms~' 
discussed in the latter half of the paper. 



From many different density functions, more or less suitable 
for modelling the distribution of Gaussian widths, the values 
of which are defined to be positive, we have chosen to use the 
lognormal distribution, as its usage was technically rather con- 
venient (if some parameter is distributed according to the log- 
normal law, the distribution of the logarithms of this parameter 
is described by a Gaussian). Mebold ( 1972) has argued that his 
narrow and shallow components are more easily distinguishable 
when the line-width distribution is studied in logarithmic scale 
and in comparison with all initially tested distributions the log- 
normal one gave the best fits. 

For the lognormal distribution, the mean value can be cal- 

culated as M(x) = ef + ~, where fi and <x may be considered 
as free parameters for fitting the model distribution to the ob- 
served one. In the case of Fig. [2^, we obtain for the means of 
three fitted distributions FWHM = 3.89 + 0.02, 11.3 + 0.1 and 
23.8 ± 0.2 kms -1 (the error estimates here are the formal errors 
of the fit and do not account for errors in determination of the 



parameters of Gaussians or any possible systematic errors due to 
the applied data reduction methods). 



The exact numerical values for the mean line-widths of the 
near zero velocity density enhancements in Fig. [2^ also depend 
on the used velocity limits (solid vertical lines in Fig.|2^). These 
limits in their turn depend somewhat on the applied latitude limit 
for the separation of the panels in Fig.[2j the lower latitude limits 
tend to indicate more symmetrical velocity limits. However, for 
reasonable changes in latitude and velocity limits the changes in 
the width estimates are rather small. For example, when lower- 
ing the latitude limit to \b\ = 30° and using the velocity limit 
|Vc| S= 8 kms -1 , we obtained the following mean width esti- 
mates: FWHM = 3.78 + 0.03, 10.8 ± 0.2 and 23.9 ± 0.5 km s" 1 . 
The comparison of these numbers with the results given above 
also yields more realistic estimates of the actual reliability of the 
mean widths than the pure formal fitting errors. 
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Fig. 3. The relation between the number of Gaussians used for 
the decomposition of the profile, the peak brightness tempera- 
ture (as estimated from the Gaussian decomposition) of the pro- 
file (red solid symbols, the vertical axis on the left) and the mean 
channel value in the region used for decomposition (blue dashed 
symbols, the vertical axis on the right). The errorbars represent 
the standard deviations around the mean of the corresponding 
distributions. The symbols are horizontally shifted to avoid over- 
lap. 



4. The Gaussians from the simplest profiles 

In the previous sections we used all the Gaussians obtained in 
the decomposition of all the profiles in the LAB database and 
tried to demonstrate how the distributions of their widths differ 
in different velocity or latitude regions. Finally, we concentrated 
our attention to the components with nearly zero velocity and 
to the relatively simple profiles at high galactic latitudes. Now 
we try to use the obtained information and limit ourselves to the 
simplest profiles and to the Gaussians describing only the low- 
velocity (relatively nearby) gas. Therefore, we must apply some 
selection criteria similar to those discussed in the previous sec- 
tion. First of all, as demonstrated by Fig. [2] we must exclude the 
complex profiles near the galactic plane. However, due to differ- 
ent complexity of the profiles in different longitude intervals the 
usage of a fixed latitude limit is probably not the best solution for 
this problem. Therefore, we decided to estimate the complexity 
of each profile individually. 

The complexity of the profile has two aspects. On the one 
hand, it is determined by the line formation process and to con- 
sider this, we must know the properties of the gas along the line 
of sight. In general, this information is not available and there- 
fore we cannot directly distinguish the lines representing the su- 
perposition of the emission of many gas concentrations projected 
on each other from those formed by the emission of a relatively 
small number of gas clouds. However, on average, the brighter 
profiles are also more complex. 

On the other hand, it is also important to know how well 
determined the parameters of the Gaussians obtained in the de- 
composition of every given profile are. In principle, this could 
be estimated. We have determined for every profile a set of 
Gaussian parameters, which corresponds to some minimal value 
of the residual rms of the deviations between the observed pro- 
file and its Gaussian representation. If we change the values of 
Gaussian parameters, the rms increases. The region within which 
the rms increases by no more than some predefined amount de- 
fines the confidence region around our solution. If we define the 



allowed increase of the rms so that the corresponding part of the 
parameter space contains some fixed percentage of the probabil- 
ity distribution of Gaussian parameters, the size and the shape 
of the confidence region give us the estimates for the errors of 
the obtained parameters. At the same time, this procedure re- 
quires huge amounts of computation, most likely not justified 
for the present study. However, on average, for profiles with a 
larger number of overlapping Gaussians the parameters of the 
latter are determined with lower confidence. Finally, as the num- 
ber of Gaussians used for decomposition is rather well corre- 
lated with the emission content of the profile (Fig. [3), we choose 
as a primary indicator of the profile complexity the number of 
Gaussians needed for its decomposition. 

However, in Paper II we demonstrated that some of the 
Gaussians are caused by different observational, reductional and 
decompositional problems and with high probability do not de- 
scribe the properties of the galactic Hi. As some profiles con- 
tain many spurious Gaussians, their inclusion in the analysis may 
considerably distort the profile complexity estimates. Therefore, 
before considering the complexity of the profiles, it is useful to 
exclude such Gaussians from consideration. We identify them by 
the criteria illustrated in Fig.[TJby thick solid lines. As described 
in Paper II, these criteria are analytically given by 



lg(FWHM) < 0.057 lg(r b0 ) 2 - 0.067 lg(T b0 ) + 0.094 



(3) 



(very narrow Gaussians are mostly due to stronger observational 
noise peaks and radio-interferences) and 



lg(r b0 ) < -0.370 lg(FWHM) 3 + 1.132 lg(FWHM) 2 - 
1.567 lg(FWHM) + 0.1 17 



(4) 



(most of the weakest components are caused by different uncer- 
tainties in the determination of the baselines for profiles). Here 
we are not interested in high-velocity clouds and therefore do not 
apply in the case of Eq. © any velocity limits used in Paper II, 
but we still reject all components with the central velocities out- 
side the velocity range used for the decomposition. After rejec- 
tion of these spurious components we classified the profiles on 
the basis of the total number of remaining Gaussians needed for 
the decomposition of any particular profile. 

So far we have used in this section all the Gaussians irrespec- 
tive of their central velocities. However, in the previous section 
we demonstrated that the mean line-widths of the components 
tend to increase with the increase of the absolute value of their 
central velocities. We also saw that this increase is slower for 
near zero velocities and becomes stronger at higher velocities. It 
is somewhat arbitrary to fix any numerical velocity limits for the 
regions with slow and faster line-width changes. In Fig.|2]these 
differences are also slightly exaggerated by the compression of 
the velocity scale at higher velocities. Nevertheless, when for the 
line-width group with widest Gaussians in Fig. [2^ in the veloc- 
ity interval of |Vc| ^ 5 kms the mean line-width changes are 
not more than 0.25 kms 1 per 1 kms 1 of the central velocity 
change, for the velocity interval of -30 ^ Vc ^ -10 km s this 
gradient increases at least to 1 .5 km s per 1 km s . Therefore, 
to restrict ourselves to the region of more or less constant mean 
line-widths, we accept for the following the velocity range ap- 
propriate for the higher latitude data and indicated in Fig. [2^ by 
solid vertical lines. After such selection we separated in each 
class of profiles the Gaussians, central velocities of which fall 
into this range of -9 < Vc < 4 km s , and constructed for every 
class the distribution of the widths of these components. 

It turned out that for profiles decomposed into one, two or 
three accepted Gaussians, the concentration of the widths around 
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Fig. 4. The distribution of the widths of Gaussians for the pro- 
files decomposed to one (upper panel), two (middle panel) or 
three (lower panel) components. 



some distinct mean values was rather well defined without any 
additional selection criteria (Fig.|4|i. To estimate the mean widths 
characteristic of these different line-width groups, we once again 
modelled the obtained distributions as sums of lognormal distri- 
butions. As in Paper II we saw that the discrimination between 
the probably real and spurious Gaussians becomes more difficult 
at larger widths, we restricted the modelling process to the line- 
widths FWHM < 100 kms . Approximately above this width 
limit the selection criterion given by Eq. (0|i becomes inopera- 



tive and at the greater widths no Gaussians are excluded by it. 
As a result, in these regions the proportion of the spurious com- 
ponents among the Gaussians, accepted for our analysis, may in- 
crease. This increase may be visible in the FWHM > 100 km s 
region of Fig. [4}: as a group of points above the fitted model line. 

As can be seen from Fig.Hk, when the profile is successfully 
decomposed by only one low velocity Gaussian, such compo- 
nents belong exclusively to the line- width group 3, as defined by 
Fig. Sj;. When the profile is decomposed by two Gaussians, the 
widths of corresponding components have already a much wider 
spread (Fig. |4j}). Most of them still fall to group 3, but also to 
group 2, first discussed in connection with Fig.|2^, is well repre- 
sented. Moreover, this is the case where group 2 is most clearly 
visible as a separate concentration of the line-widths. Beside 
these two dominating groups two weaker ones can be seen. On 
the left-hand side of the distribution, signs of group 1 , which was 
the main group in Fig.[2j3, become visible, but on the right-hand 
side the presence of a new group seems to be rather obvious. 

In Fig. H}; the role of group 1 has considerably increased 
and it has become a well-established component of the distri- 
bution. The new fourth group still remains rather ill-defined. 
Moreover, as with more complex profiles the role of extremely 
wide Gaussians increases, the presence of group 4 becomes even 
somewhat more questionable than it may seem on the basis of 
Fig- Eh. Therefore, we add this new group to our list of the pos- 
sible frequency enhancements in the distribution of the Gaussian 
widths, but we must remember that this group is established con- 
siderably less reliably than the other three. 

On the basis of the profiles decomposed into one, two or 
three low velocity Gaussians, the estimates for the mean line- 
widths of the four groups of the components described above, are 
4.9 ±0.1, 12.0 + 0.3, 24.4 + 0.2, and 43 ±5 kms" 1 . Unlike the pre- 
vious estimates, here the errors are not given on the basis of the 
formal errors of the parameter fitting, but a procedure somewhat 
similar to bootstrap was used. We estimated the mean values and 
their formal errors separately for profiles with one, two or three 
Gaussians and also for any combination of such classes (the pro- 
files with one and two or one and three Gaussians together and 
so on), calculated the weighted (on the basis of formal fitting er- 
rors) means of the results and accepted as an error the standard 
deviation of the obtained estimates from these means. We expect 
that such estimates include more realistically the uncertainties of 
the obtained values than the pure formal errors of the fit of the 
model to a fixed dataset. 

For more complex profiles (four and more accepted 
Gaussians per profile) the components of the first group (the nar- 
rowest ones) become more and more dominating and this makes 
the recognition of the second group harder and harder (this is 
well demonstrated also by Fig. [2] where the second group is in- 
visible in the lower panel, which represents the complex pro- 
files near the galactic plane). The third group remains relatively 
well established, but the identification of the fourth group be- 
comes increasingly complex and dependent on the accepted up- 
per width limit, as in such profiles the frequency of extremely 
wide components increases. However, these changes are weaker 
for weaker profiles. At the same time, there is no clear-cut line 
between weak and strong profiles. Therefore, we decided to at- 
tempt to determine the width distributions of the Gaussians for 
different profile strength limits and we used in the following two 
different ways to define the profile strength. In the first case, 
we used as an indicator of the profile strength the peak bright- 
ness temperature of the profile (as estimated from the smooth 
Gaussian representation), but in the second case, we used the 
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Fig. 5. The sky positions of the 79 648 profiles 
used for the final line-width estimates. The po- 
sitions of the simplest profiles are marked by 
the reddest dots. The colour scale represents in 
the first order the number of Gaussians per pro- 
file. Among the profiles with an equal number 
of Gaussians those with lower peak brightness 
temperature are plotted with redder colour. 



mean channel value (calculated from the original observed pro- 
file) over the velocity range used for decomposition. 

For both profile strength definitions we repeated the esti- 
mates over a wide range of strength upper limits, but accepted 
for final results only the cases, where the means for all lognormal 
distributions were estimated with the formal errors below 50% 
(or 25% in a more restrictive case). In the case of the strength 
definition through the peak brightness temperature, the profile 
strength upper limit was increased in steps of 0.1 K from the 
zero until there were no more acceptable results during the last 
10 steps. Only for the profiles decomposed into three Gaussians 
even the inclusion of the strongest lines with the highest peaks 
extending up to 44.3 K permitted determination of the parame- 
ters of different line-width groups. In all other cases no profiles 
with Tb.max > 15.8 K could be used and 89.7% of the profiles 
used for the final estimates had Tb.max < 12 K. 

When using the peak brightness temperature as the profile 
strength indicator, the obtained mean line-widths were 4.3 ± 0.6, 
12.0 ± 0.5, 24.1 ± 0.7, and 42 ± 5 km s" 1 for the 50% maximal 
error limit and 4.6 ± 0.5, 12.1 ± 0.4, 24. 1 ± 0.6, and 41 ± 4 km s _I 
for 25%. Then the same procedure, but increasing the profile 
strength upper limit in steps of 0.001 K was also repeated for the 
mean channel value as a profile strength indicator. In this case, 
the results were 3 .9 + 0.6, 1 1 .7 ± 0.4, 24. 1 ± 0.6, and 43 ± 5 km s~ 1 
for the 50% limit and 4.3 ± 0.7, 1 1 .8 + 0.4, 24. 1 + 0.7, and 42 ± 
4kms"' for 25%. 

As we can see, the results are in general agreement within 
their error estimates, which consider the dispersion of the results 
for a particular selection process, but do not take into account 
the uncertainties related to the choice of the peak brightness or 
mean channel value as the profile strength indicator or the usage 
of the 25% or the 50% precision limit. As the results seem to de- 
pend more on the way of estimating the profile strength than on 
the limit of formal errors of the parameter determination, we ac- 
cept as the final result the average over different profile strength 
estimates for the 50% error limit. After correcting the mean line- 
widths also for the width of the correlator channels in the LAB 
the corresponding values are: 3.9 + 0.6, 11.8 + 0.5, 24.1 + 0.6, 
and 42 ± 5 kms -1 . Altogether these estimates are based on the 
decomposition of 79 648 selected profiles. This set of profiles 
represents 42.0% of all observed sky positions in the LAB. The 
sky distribution of the corresponding profiles is given in Fig. [5] 



As an additional check, we repeated the procedure described 
in this section also for the symmetrical velocity limit (|Vcl - 
8 km s _1 ) of the analysed profiles. As a profile strength indica- 
tor, we used the peak brightness and accepted the fits where the 
means for all lognormal distributions were estimated with the 
formal errors below 50%. In this case, we obtained the follow- 
ing values: 4.3 ± 0.5, 12.0 ± 0.5, 24.0 ± 0.7, and 43 ± 5 kms" 1 , 
the results still in good agreement with the earlier ones. 

As we can see from Fig.|5j when selected on the basis of the 
number of Gaussians per profile and the strength of the emis- 
sion line, the profiles used do not follow any well-defined lati- 
tude limit. Nevertheless, the mean values obtained for different 
line-width groups are still in good agreement with the estimates 
from the latitude limited samples obtained in Sec. 3. (we remind 
that the error estimates in these two cases are based on different 
considerations). This result, when combined also with Fig. [2] 
indicates that when there are clear differences in the overall line- 
width distribution of the most simple and most complex profiles, 
the transition from one to another is relatively smooth, without 
a well-defined boundary between them. Moreover, as differently 
selected samples give rather similar values for the mean line- 
widths of the groups, it seems that the width distributions of the 
profiles of different complexity differ from each other mainly by 
the degree of the presence of different line-width groups and not 
so much by the mean widths of these groups. 

5. The comparison with earlier results 

In the previous section we demonstrated that in the distribution 
of the widths of the Gaussians representing the simplest 21 -cm 
line profiles of the galactic H i, we can distinguish at least three, 
but probably even four distinct groupings of the line-widths. 
Earlier, similar results have been reported by Verschuur and co- 
authors, who have published several papers (119891 119941 119991 
hereafter VP, [20021 [20041 hereafter V4) where they argue that 
neutral hydrogen emission profiles produced by gas in the lo- 
cal interstellar medium are characterised by three, or probably 
four, line-width regimes where dominant and pervasive features 
have widths of the order 5.2 + 0.4 (VP regime 3), 12.8 + 0.4 (VP 
regime 2), 31.0 + 0.5 (VP regime lb), and 50. 1 ± 0.6 km s" 1 (VP 
regime la). They also note that these line-width regimes show 
a striking resemblance to a set of velocity regimes described by 
a plasma physical mechanism called the critical ionization phe- 
nomenon. 
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As we can see, our results are qualitatively identical. We 
both clearly distinguish three different groupings or regimes of 
the H i line-widths and admit that there may be also the fourth 
one with a larger mean line-width than the three main groups. 
However, we must recognize that there are considerable differ- 
ences in the numerical values of the mean line-widths obtained 
by us and by Verschuur (as given in VP). On average, our re- 
sults are about 17% smaller than those given in VP. The dif- 
ferences are the largest (approximately 22%) for groups 1 and 
3 (VP regimes 3 and lb) and the smallest (12%) for groups 2 
and 4 (VP regimes 2 and la). At the same time, there are also 
considerable differences in our approaches to data analysis. Our 
studies are based on somewhat different initial data (the LAB 
instead of the published version of the LDS), we have used dif- 
ferent Gaussian decomposition algorithms and selection criteria 
for choosing the components for final analysis. Our approaches 
to estimating the mean line-widths of the groups of the compo- 
nents are also different. Therefore, it is interesting to check how 
all these differences may have influenced the results obtained. As 
in most cases in his papers Verschuur has not given the results 
for his regime la, we follow the same practice in this section. 

In Gaussian decomposition Verschuur and co-authors have 
clearly preferred the human assisted approach, as they assume 
that Gaussian decomposition "is not something that can be left 
to an idealized computer program" (VP). We have relied on 
the computer from the beginning (Paper I). The probability that 
the differences in decomposition algorithms may be the primary 
source of the discrepancies in the mean line-widths of our line- 
width groups and Verschuur's regimes could be reduced if we 
could demonstrate that our programme gives for identical ini- 
tial data practically the same results as obtained by Verschuur. 
This can be most easily done for the "random directions cho- 
sen from the Leiden/Dwingeloo data files" for which VP give in 
their Table 2 the list of sky positions from where the profiles of 
the sample have been chosen. 

We analysed the Gaussians corresponding to these profiles 
of the LDS. In the width range of 0.5 < FWHM < 60 kms" 1 , 
used by VP (see Fig. 4b in VP), we found three concentrations 
of the line-widths around mean values of FWHM — 5.4 + 0.7, 
13.0 ± 0.5 and 29.7 + 6.8 kms" 1 . In their Table 7, VP gave for 
the same sample the values equal to 5.2 ± 1.3, 12.8 ± 1.2 and 
32.8 ± 1.8 kms" 1 . As can be seen, for exactly the same initial 
data our results coincide within errors (here we use once again 
the formal errors of the fit). 

The general selection criteria used by Verschuur are best de- 
scribed in the second chapter of V4. To test these criteria and 
to approximate the V4 results, we decomposed 206671 profiles 
of the LDS (as published by Hartmann & Burton 1997) into 
1 644 665 Gaussians and applied the following selection crite- 
ria defined in accordance with our understanding of the criteria 
used in V4: 

1. Profiles along the lines of constant Galactic longitude at 10° 
intervals, starting with I = 0° and every 1 in latitude were 
used. 

2. The latitude range was from \b\ = 80° (inc.) down toward the 
Galactic plane until: 

(a) the latitudes \b\ < 40° were reached, or 

(b) the unsmoothed profile was found in which the bright- 
ness temperature in at least one channel in the velocity 
range of -460 < Vc < 396 km s was 10 K or greater. 

3. If the smooth profile, as reconstructed from the Gaussian 
decomposition, contained in the velocity interval of -60 < 
Vc < 40 kms -1 multiple peaks with separation less than 



20 km s and the depth of the minimum between them less 
than 10% of the height of the lower peak, the Gaussians cor- 
responding to this profile were not used in the analysis. 

4. Any Gaussian with a peak brightness temperature of < 
0.15 K was removed from further consideration. 

5. Only the components with the central velocity Vc in the 
range of -31 < Vc < 18 km s were used for the study and 
all profiles requiring more than four such Gaussians were re- 
moved from the analysis. 

Following V4, we analysed the results separately for positive 
and negative latitude data. For positive latitudes we obtained the 
mean values FWHM = 4.5±0.5, 14.8+1 .7, 30.4+ 1 .9 km s" 1 and 
for negative-latitudes we found FWHM = 3.9 + 0.5, 1 1.5 ± 1.9, 
23.2 + 0.4 kms _I . We can see that on average these numbers 
are slightly larger than the corresponding results based on our 
selection criteria and the LAB data. Most likely this is due to 
the wider velocity range used by V4 (indicated in Fig. [2^ by two 
dashed vertical lines). However, in the first order, these numbers 
must be compared with those given in Table 7 of V4: FWHM = 
5.8 + 2.4, 13.7±2.3,35.1±4.0kms _1 for positive and FWHM = 
5.2 ± 1 .9, 14.5 ± 2.4, 33.9 + 3.6 km s -1 for negative latitude data, 
respectively (here our and V4 error estimates are most likely not 
directly comparable). With the exception of group 2, our values 
are still smaller than those by V4. 

So far we have estimated the mean line-widths of the differ- 
ent line-width groups through the approximation of the distribu- 
tion by the sum of lognormal functions. At the same time, the 
mean values given in Table 7 of V4 are calculated as averages 
over the predefined ranges of widths. If we apply the same ap- 
proach and the same line-width ranges to our results described in 
this chapter, our mean line-widths are for positive-latitude data 
FWHM = 5.7+2.2, 14.3±2.8,34.4±4.7kms _1 and for negative- 
latitude data FWHM = 5.6 + 2.2, 14.2 + 3.0, 34.4 ± 4.9 kms" 1 
(here the standard deviations of the line-widths around the means 
in predefined width intervals are given as error estimates as 
seems to have been done also by V4). These numbers are in good 
agreement with the V4 results as cited above. However, in this 
case it remains somewhat questionable how justified the averag- 
ing ranges used in V4 are. For example, why are in Table 7 of 
V4 the Gaussians with widths of 20 < FWHM < 28 km s _1 not 
included into any line-width regime? 

6. Discussion 

In V4 the obtained mean line-widths are used to estimate the 
kinetic temperatures of the gas. Verschuur found that the tem- 
peratures corresponding to the line-width regimes lb and 3 are 
24 000 K and 750 K, respectively. These temperatures were con- 
sidered as physically unreal. On the basis of the coincidence of 
the obtained mean line-widths with the critical ionization veloc- 
ities (CIV) for some atomic species (6 kms -1 for Na and Ca, 
13.5 kms" 1 for C, N, and O, 34 kms" 1 for He and 51 kms" 1 
for H, respectively) a model with collisionless gas, in which the 
origin of the profile broadening represents the influence of the 
critical ionization velocity effect, is presented. This interpreta- 
tion has not found wide acceptance in astronomical community 
and our results make such explanation even more unlikely as the 
agreement between the selected CIVs and the mean line-widths 
seems to be poorer than stated by VP. 

However, some decades ago, Clark d 19651 1 proposed a 
two-component model for interstellar Hi where cold absorb- 
ing clouds are surrounded by hot, ubiquitous intercloud gas. 
Theoretical studies (e.g. Field et al. 1 19691 1 also demonstrated that 
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the local H i could be considered as a two-phase medium, where 
much of the gas is either a warm neutral medium (WNM) with 
T ~ 8000 K or a c old ne utral medium (CNM ) with T ~ 50 K 
( Kulkarn i & Hei les [T9871 Dickey & Lockman [19901 Wolfire et 
al. |1995i Ferriere 1200 11 1. At present this is one of the most widely 
accepted models for the local H i. Therefore, we try to check how 
our four line-width groups may be related to this model. 

Mebold ( |19721 l has estimated that the emission line compo- 
nents, corresponding to the WNM, have a Gaussian dispersion 
of about 9.7 km s or FWHM w 23 km s _1 . This width is rather 
close to the mean width of our third line-width group and there- 
fore we may expect that this line-width group may be related 
to the presence of the WNM in the ISM. With other line-width 
groups the situation seems to be somewhat more complicated. 

Verschuur ( 119741 1 and Schwarz & van Woerden ( 119741 1 have 
demonstrated that the histogram of the CNM H i emission line 
widths peak around <r v = 2.2 kms~', which corresponds to 
FWHM 5.2 kmT 1 - a value somewhere between our first 
and second line-width groups and directly not comparable to any 
line-width group. However, here we must notice that their results 
have been obtained using all narrow emission line components 
presented in the observed profiles. At the same time, as discussed 
in the Introduction, for stronger emission components the satu- 
ration effects become important and this makes the line shape 
nongaussian and such lines seem to be broader than unsaturated 
profiles. 

In our analysis we have used only relatively weak lines with 
nearly 90% of the used profiles having peak brightness less than 
12 K. By combining the equation of radiative transfer with the 
approximate correlation between the H i spin temperature and 
optical depth t (e.g. Payne et al. 1983), we may estimate that for 
most of our profiles r < 0.1, and therefore, we may use the ap- 
proximation for optically thin gas T\,(v) = Tst(v). For r < 0.1 
this approximation introduces up to 5% error in the peak bright- 
ness temperature of the emission profile, but practically does not 
change the line-width and consequently we must compare our 
line-widths not with those of much stronger emission lines cited 
above, but with the line-widths of the absorption lines of the 
CNM. 

Crovisier ( 1981 ) has reported that the distribution of the ve- 
locity dispersion inside the CNM clouds peaks around 1 kms~', 
but the mean value for the distribution is about 1.7 kms -1 . The 
latter corresponds to FWHM « 4.0 kms~', which coincides 
within errors with the mean line-width of our first line-width 
group. This makes it rather plausible that the existence of both 
line-width groups 1 and 3 is a reflection of the presence of two 
thermal phases of neutral gas in the local interstellar medium. 

Theoretical CNM and WNM temperatures are derived by 
calculating the equilibrium temperature as a function of thermal 
pressure. As pointed out above, there exist two stable ranges of 
equilibrium, separated by a region of unstable temperatures. If 
we interpret the line-width of group 2, 11.8 kms~', as thermal, 
the corresponding kinetic temperature would be about 3000 K, 
which lies just in the region of unstable temperatures. At the 
same time, the presence of some gas with such temperatures 
seems to have been a long-standing problem. Cox (2005) has 
recently given a review of the three-phase model of the interstel- 
lar medium, but there the additional third phase (warm H ia) has 
a temperature of 5000 K. Heiles (2001 ) and Heiles & Troland 
(120031 1 have reported that at least 48% of the WNM lies in the 
thermally unstable region of 500 - 5000 K. Similar results were 
obtained by Mebold et al. ([1982 ). On the basis of the results 
by Dickey et al. JT977l [19781 [19791 , Crovisier ( fl98lT > discussed 
the presence of cold clouds, lukewarm clouds and not strongly 



absorbing material. Therefore, with some probability our line- 
width group 2 may also represent thermally unstable neutral gas 
and our (also Verschuur's) results (Fig.|2^ and especially Fig.[4}3) 
demonstrate that this is not just a tail of the CNM or the WNM 
temperature distribution, but a component in its own rights. 

Besides first three rather obvious line-width groups there is 
also the line-width group 4 with the mean line-width around 
42 km s . As mentioned above, the real existence of this group 
is much more questionable, as it may be seriously affected 
by baseline and stray radiation corrections in the survey (see 
Paper 2). Therefore, we make only some brief remarks about this 
group. Kalberla et al. ([1998 ) have argued in favour of the exis- 
tence of some neutral gas in the galactic halo. The mean velocity 
dispersion of this gas is a — 60 ± 3 km s~ at the north galactic 
pole and it increases at lower latitudes. The corresponding Hi 
line-widths must be more than 3 times higher than the mean for 
our line-width group 4. Line-widths (with <r — 21 + 4 kms 1 ) 
much more similar to our group 4 are reported from the measure- 
ments of absorption lines by Dwarakanath (2004) and Mohan et 
al. (I2004all2004bl 

It is also interesting to estimate the mass fractions of the 
Hi gas, corresponding to the described four line-width groups. 
However, in our case these estimates are made rather question- 
able by at least two factors. First of all, the line components with 
highest peak brightness temperatures belong to the line-width 
group 1 and in the following groups the components are on aver- 
age progressively weaker. As for our estimates we have selected 
only the simplest and the weakest profiles, we have biased the 
mass fraction estimates in favour of the groups with wider and 
weaker line components. Therefore, for the estimates of mass 
fractions, it is more natural to return to our latitude limited sam- 
ples (similar to Fig. |2^), but even in this case the stronger lines 
near the galactic plane become optically thicker, which causes 
the underestimation of corresponding masses. 

The second problem is that when the mean line-widths of 
our line-width groups are well determined by log-normal fits, the 
standard deviations of the corresponding distributions are rather 
poorly constrained due to correlations between these parameters 
for different line-width groups. To obtain the fits for different lat- 
itude limits \b\ > £>ii m , we had therefore to restrict the variability 
of these parameters and this probably introduces some system- 
atic errors into the mass fraction estimates, specially for lower 
values of b\i m . 

Nevertheless, we found that for \b\ > 40° the velocity inter- 
val of -9 < Vc < 4 kms _I contains more than 51% of all Hi. 
As the line-width groups found in this velocity interval seem 
to continue even beyond the accepted velocity limits (with the 
exception of group 2, which is visible only in the vicinity of 
Vc = kms~', they just become smoothly wider at higher ve- 
locities), we may conclude that the described line-width concen- 
trations are characteristic of most of the high-latitude H i. 

The distribution of the H i between the obtained line-width 
groups depends on the latitude limit used. Only the relative con- 
tent of group 4 is about 4% for all samples \b\ > bu m with dif- 
ferent values of b\\ m . For 15° <: b\\ m <: 55° the relative fractions 
of other groups change rather slowly and smoothly. The share of 
groups 1 and 2 decreases from about 18% and 45% for \b\ > 15° 
to about 11% and 39% at \b\ > 55°, respectively. However, we 
must stress once more that the estimates for these groups are 
strongly correlated and especially for lower galactic latitudes the 
actual fraction of group 1 may be higher and that of group 2 
lower than stated here. The content of group 3 increases for the 
same latitude range from about 34% to 46%. For larger values 
of b\i m the relative fraction of group 1 remains at about 11%, but 
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the share of group 2 quickly increases to about 47% at bu m = 75° 
and that of group 3 drops to about 37%. Therefore, for most of 
the relatively high-latitude gas we, like Heiles (2001), may con- 
clude that about half of the WNM is represented by line-width 
group 2 which temperatures lie in the thermally unstable regime. 
For b]im < 15° our analysis is most likely not applicable at all 
and for by im > 75° the results become uncertain due to the small 
number of profiles. 

7. Conclusions 

We have demonstrated that for simple (Sec. 4.) and high-latitude 
(Sec. 3.) profiles the line-widths of the Hi emission profiles ex- 
hibit concentration to certain preferred line-width groups. With 
this we confirm and extend the earlier similar results obtained 
by Verschuur and co-aufhors ( [19891 [19941 [19991 [20041 . We de- 
fine the limits of the usable velocities (-9 < Vc < 4 kms~') 
on the basis of the extent of the region, where the dependence 
of the line-widths on the component velocities is the weakest. 
We have not fixed any exact criteria for the simplicity of the 
profiles, but have tried instead many different combinations of 
the limits on the number of the Gaussians in the decomposition 
and on the profile strength as defined by their peak brightness 
or mean channel value. For the final results we have accepted 
only these combinations of the limits, which allow the determi- 
nation of the mean line-widths of all width groups with the for- 
mal fitting errors less than 50%. This condition was fulfilled for 
all profiles, decomposed into one, two or three components. For 
such profiles we demonstrated that it is possible to distinguish 
three or four groups of preferred line-widths. The existence of 
groups 2 and 3 is very clearly visible, when using only the pro- 
files decomposable by one or two Gaussians (Figs. [4^ and[4j}). 
Group 1 can be better distinguished in profiles decomposed by 
three components (Fig. [4}:). The reality of group 4 remains more 
questionable. A similar separation was possible also for weaker 
profiles decomposed into four, or in some cases even into up 
to 7 Gaussians. The mean line-widths of the groups found are 
3.9 + 0.6, 11.8 ±0.5, 24.1 ±0.6, and 42 ± 5 km s" 1 . 

The mean line- widths obtained by us are somewhat smaller 
than those proposed earlier by Verschuur. The disagreement may 
be caused partly by different algorithms used to measure the 
mean line-widths of our groups and VP regimes. Another im- 
portant point appears to be the usage by V4 of Gaussians from a 
considerably wider velocity range than it is acceptable to us. Due 
to the smaller mean line-widths of our line-width groups, also 
the interpretation of the existence of the corresponding groups 
through the collisionless gas and critical ionization velocities, as 
proposed by Verschuur, seems rather questionable to us. More 
likely groups 1 and 3 represent CNM and WNM components of 
the ISM. The existence of group 2, which may represent up to 
40% of all Hi, in the thermally unstable regime, is in conflict 
with usual two component equilibrium models of the ISM, but 
from our data (Figs. [2^ and [4}}) the existence of such a group 
seems to be well established. The reality and the origin of the 
broad lines with the widths of about 42 km s is more obscure. 
These, however, contain only about 4% of the total observed col- 
umn densities. 
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